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Abstract 

The  spherical-grain  mathematical  model  is  modified  and  tested  against  experimental  data  for  single-layer,  gas-diffusion 
electrodes  of  alkaline  fuel  cells.  The  model  assumes  that  the  electrode  is  made  of  spherical  agglomerates  of  Raney  metal 
and  polytetrafluoroethylene  (PTFE)  that  are  flooded  with  electrolyte;  the  gas  occupies  the  macropores  of  the  electrode.  In 
addition  to  previous  analysis  of  the  diffusion  and  reaction  in  the  grains,  the  modified  model  includes  the  resistance  of  gas 
diffusion  into  the  macropores  and  a  thin  electrolyte  film  surrounding  the  grain.  The  original  model  and  the  modified  model 
are  both  compared  with  experimental  polarization  data  for  hydrogen  oxidation  on  an  Ni/PTFE  electrode  in  alkaline  electrolyte. 
The  newly  developed  model  predicts  accurately  the  experimental  data  in  all  regions. 
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1.  Introduction 

In  the  authors’  laboratory,  a  dry  method  is  used  to 
produce  porous  gas-diffusion  electrodes  for  fuel  cells. 
These  electrodes  have  only  one  layer,  unlike  most  fuel¬ 
cell  electrodes  that  have  two  layers,  namely,  a  gas- 
diffusion  and  an  active  layer.  Numerous  attempts  have 
been  made  to  model  fuel-cell  electrodes  mathematically 
in  order  to  obtain  theoretical  polarization  behaviour 
and,  thus,  to  optimize  their  structure  and  composition. 
Conventional  cylindrical  models,  if  applied  to  the  single¬ 
layered  electrodes,  will  not  show  the  limiting-current 
behaviour  due  to  the  absence  of  a  gas-diffusion  layer 
and,  thus,  to  the  absence  of  a  gas-diffusion  resistance 
[1],  A  thorough  review  regarding  mathematical  models 
for  electrode  polarization  has  been  given  by  Celiker 
et  al.  [2,3].  In  these  publications,  a  new  model  has 
been  proposed.  This  is  called  ‘the  spherical-grain  model’ 
in  which  the  electrode  is  considered  to  have  a  single 
layer.  The  model  assumes  (based  on  surface  analysis) 
that  Raney  metal  particles  are  dispersed  in  the  electrode 
as  spherical  grain  that  are  flooded  with  the  electrolyte, 
while  the  macropores  (i.e.,  the  space  between  the  grains) 
are  occupied  by  the  gas  phase.  The  spherical-grain 
model  predicts  that  the  current  density  will  be  almost 
linear  with  the  overpotential.  The  previous  model  has 
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been  compared  with  experimental  data  obtained  by  Al- 
Saleh  et  al.  [4]  for  a  Raney  Ni/PTFE  electrode.  The 
model  agrees  with  the  experimental  data  only  at  low 
temperatures  and  low  overpotentials.  At  higher  over¬ 
potentials  and  temperatures,  part  of  the  experimental 
data  deviates  from  the  model  and  indicates  the  approach 
to  a  limiting  current.  This  discrepancy  can  be  explained 
on  the  basis  of  neglected  external  mass-transfer  resis¬ 
tances.  In  this  paper,  the  spherical-grain  model  is 
modified  by  incorporating  the  mass-transfer  resistance 
in  the  macropores  for  gas  diffusion  and  by  assuming 
the  existence  of  a  thin  film  of  electrolyte  around  the 
spherical  grains. 


2.  Modified  flooded  spherical  agglomerate  model 

The  electrode  is  assumed  to  be  made  of  porous, 
spherical,  catalyst  grains  that  are  held  together  by  fine 
PTFE  strands.  The  space  between  these  spheres  con¬ 
stitutes  the  macroporosity.  The  presence  of  the  PTFE 
strands  renders  these  macropores  hydrophobic.  There¬ 
fore,  no  electrolyte  remains  in  the  macropores.  The 
gas  has  first  to  diffuse  into  the  macropores  before  it 
reaches  the  grain.  In  this  modified  model,  the  resistance 
associated  with  macropore  diffusion  has  been  included. 
Besides  this,  it  is  proposed  that  the  hydrophilic  grain 
is  surrounded  by  a  stagnant  film  of  electrolyte  of  constant 
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thickness.  The  reactant  gas  is  dissolved  at  the  outer 
surface,  and  then  diffuses  through  the  stagnant  film  to 
reach  the  surface  of  the  flooded  catalyst  grain.  Further 
diffusion  of  these  molecules  takes  place  inside  the  grain 
with  simultaneous  reaction.  The  ionic  conductivity  path 
(explained  previously  [3]  in  terms  of  narrow  electrolyte 
bridges  between  the  grains)  is  more  convincing  now 
since  the  electrolyte  film  is  taken  to  surround  the  grain. 
The  working  mechanism  and  the  structure  of  the  elec¬ 
trode  are  shown  schematically  in  Fig.  1.  The  resistances 
in  the  different  sections  of  the  diffusion  path  of  the 
gas  are  also  shown. 

For  the  mathematical  development  of  the  modified 
model,  it  is  necessary  to  make  the  following  assumptions: 

(i)  Raney  catalyst-PTFE  grains  are  spherical,  iden¬ 
tical  in  size/structure  activity,  isopotential,  distributed 
homogeneously  in  the  electrode,  flooded  with  electro¬ 
lyte,  and  touching  one  another. 

(ii)  There  exists  a  thin  layer  of  the  electrolyte  on 
the  surface  of  each  agglomerate.  This  film  helps  in  the 
ionic  and  electronic  transport  across  the  dissolved  gas. 

(iii)  The  local  current  density  in  the  grain  is  given 
by  the  Butler-Volmer  equation. 

(iv)  Diffusion  of  the  dissolved  gas  in  the  grain  and 
the  film  occurs  in  the  radial  direction  only,  while  the 
macropores  diffusion  is  in  an  axial  direction  only. 


(v)  The  potential  changes  only  across  the  thickness 
of  the  electrode. 

(vi)  There  is  no  kinetic  limitation  on  the  dissolution 
of  gas  in  the  electrolyte. 

(vii)  The  ionic  concentration  equilibration  along  the 
electrode  is  achieved  by  an  evaporation/condensation 
mechanism. 

The  modified  model  is  based  on  analysis  of  the  grain, 
the  film  and  the  macropore.  These  analyses  are  then 
used  with  Ohm’s  law  to  describe  the  relationship  be¬ 
tween  current  density  and  overpotential. 

2.1.  The  grain 

In  the  grain,  the  dissolved  gases  are  diffusing  radially 
and  reacting  simultaneously.  The  local  current  density 
in  a  grain  at  plane  x  in  the  electrode  is  given  by  the 
Butler-Volmer  equation: 

.  .[C  (azFV(x)\  ( ~(l-a)zFV(x)\l  _ 

,-‘°[c0expl-vrTM,pl — — )\  (1> 


The  rate  of  consumption  of  the  reactant  gas  in  the 
grain  in  terms  of  mol  cm-3  s-1  is  given  by: 


Fig.  1.  Schematic  diagram  of  fuel-cell  electrode  and  Raney  Ni/PTFE 


The  mass  conservation  equation  with  simultaneous  dif¬ 
fusion  and  reaction  under  steady-state  conditions  is 
given  by: 


Combining  Eqs.  (1)  to  (3)  with  necessary  rearrangements 
yields: 


where, 


A-i  =  exp 


azF-nix) 

RgT 


Aj  —  exp 


~(l-a)zFv(x) 

RgT 


The  boundary  conditions  are: 
r=R  C  =  CS 


The  analytic  solution  of  Eq.  (4)  is: 


(5) 

(6) 
(7) 
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C(x,  r)  \CS  A3 1  R  smh(rMyA2/C~)  ^ 
C0  ~[C0~  A2\  r  sinh (RslAxA2IC0)  +  2l2 


2.2  The  film 


With  the  presence  of  an  electrolyte  film  around  the 
spherical  grains,  the  dissolved  gas  must  diffuse  first 
into  the  film  surrounding  the  spherical  grains.  If  the 
concentration  at  the  outer  surface  of  the  film  is  kp(x), 
while  at  the  interface  of  film  and  grain,  is  Cs,  then 
the  continuity  equation  for  the  diffusing  system  at  steady 
state  is  given  by: 


with  the  boundary  conditions: 

r=R  Cf=Cs 

r^R  +  8  Cf=kp(x) 

The  solution  of  Eq.  (10)  with  these  boundary  conditions 

is: 


'3Cf\  _  R(kp(x)  -  Cs) 

tfr/r_*+8  8(R  +  S ) 


(11) 


At  the  interface  between  the  film  and  the  grain,  the 
following  mass  balance  holds: 

4ir(/?  +  5)2A^fj  =477/?2Z)le^j  (12) 

Substituting  and  from  Eqs.  (9) 

and  (11)  into  Eq.  (12)  yields: 


d 2£ 


(14) 


where  Dg  is  the  effective  gas  diffusivity  in  the  macropores. 
The  term  R*  is  the  rate  of  consumption  per  unit  volume 
of  the  electrode  and  is  given  by: 


,  (*£\  IK' 

lelar/JL. 


where, 

(i-w) 
ttR2 


Combination  of  Eqs.  (14)  to  (16)  gives: 

=  0 


dp  _  4£>,e(l  - e)  /0C\ 
Ds  Hr2  I.  \  ar  I 


(15) 


(16) 


(17) 


with  the  following  boundary  conditions: 


x=0 

x=L 


dr 

P=Po 


2.4.  Ionic  overpotential 


The  overpotential  inside  the  electrode  varies  in  the 
axial  direction.  The  current  density  and  overpotential 
relationship  are  related  by  Ohm’s  law,  namely: 


d27  =  -;(*) 

dr  vR2KeN1 


(18) 


where  N j  is  given  by  Eq.  (16).  The  current  density 
generated  in  a  differential  axial  distance  dr  is: 


-Aj{x)  =  2trRN,nFDu 


dr 


(19) 


Differentiating  Eq.  (18)  gives: 
d2r?  _  1  ( dj(r)\ 

dr2  Jte(l-e)\dr/ 

Substitution  of  from  Eq.  (19)  into  Eq. 
yields: 

d2T?  =  2 nFDle  /  dC) 
dr2  kcR  \  dr  j  _R 


(20) 

(20) 

(21) 


2.3.  The  macropores 

In  the  hydrophobic  macropores,  the  reactant  gas  is 
diffusing  and  is  being  consumed  simultaneously  under 
the  following  mass-conservation  conditions: 


with  the  following  boundary  conditions: 
x=0  V  =  Vo 
d V 


x=L 


dr 


=  0 
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Substituting 


into  Eq.  (21)  produces: 


d2V  2nFD]eC0[Cs  Af 
dr2  ~  kcR  [c0  A2 


To  obtain  polarization  data  for  the  electrode,  Eqs. 
(17)  and  (22),  with  their  boundary  conditions,  have  to 
be  solved  simultaneously.  An  analytic  solution  of  these 
equations  is  not  possible.  Therefore,  a  numerical  tech¬ 
nique  is  used.  In  the  latter,  the  given  boundary-value 
problem  is  converted  into  an  initial-value  problem  (IVP) 
by  guessing  the  initial  conditions  that  are  not  given. 
Other  sets  of  initial-value  problems  are  generated  by 
differentiating  the  first  set  with  respect  to  the  assumed 
values  of  initial  conditions  [5].  All  these  sets  are  solved 
simultaneously  using  the  Runge-Kutta  method  with 
adjustable  step  size  in  double  precision. 


3.  Results  and  discussion 

The  developed  model  was  tested  in  two  parts.  In 
the  first  part,  a  base  set  of  the  parameters  was  selected 
by  using  hydrogen  oxidation  on  Ni/PTFE  in  25%  KOH 
electrolyte.  Out  of  these,  only  one  parameter  at  a  time 
was  varied  and  polarization  curves  were  plotted.  The 
observed  effect  of  that  particular  parameter  was  com¬ 
pared  with  the  theoretically  expected  trend.  Individual 
studies  were  made  of  exchange-current  density,  grain 
radius,  macropore  diffusivity,  diffusivity  in  the  electrolyte 
and  film  thickness.  In  the  second  part  of  the  test,  the 
model  equations  were  solved  for  different  temperatures 
and  plotted  against  experimental  polarization  data  for 
hydrogen  oxidation  on  Ni/PTFE  electrode  in  25%  KOH 
electrolyte.  The  same  experimental  data  were  also 
compared  with  the  previously  developed  flooded  spher¬ 
ical-grain  model. 

3.1.  Effect  of  parameters  on  polarization 

The  prediction  of  the  modified  spherical-grain  model 
has  been  investigated  by  considering  the  anodic  oxi¬ 
dation  of  hydrogen  in  alkaline  electrolyte  on  an  Ni/ 
PTFE  electrode: 

H2(aq)  +  20H"  (aq)  ^  2H20(1)  +  2e"  (23) 

The  effect  of  exchange-current  density,  diffusivity  of 
hydrogen  in  both  macropores  and  electrolyte,  grain 
size,  and  film  thickness  was  examined  by  changing  one 
parameter  at  a  time  and  fixing  the  others.  For  this 
purpose,  the  estimated  and  fitted  parameters  at  75  °C 
from  Table  1  were  selected  to  test  the  prediction  of 


the  modified  model  for  the  performance  of  a  Ni/PTFE 
electrode  in  25%  KOH  electrolyte.  The  parameters 
were  as  follows: 

Ke  =  0.464  ft1  cm  '1 

k  =  0.507  X  10"9  mol  cm-3  kPa"1 

Dg=(D'g)fc  =  2.279XlO"11  mol  cm"1  kPa"1  s"1 

£>,c  =  9.529x10" 5  cm2  s"1 

S=0.8  pm 

i0  =  1  pA  cm"2 

p0=138  kPa 

The  effect  of  the  film  thickness  (6)  was  investigated 
by  varying  it  in  the  range  from  0.05  to  1  pm.  It  can 
be  observed  from  Fig.  2  that,  with  low  values  of  8,  the 
current  generated  for  a  given  overpotential  is  high  due 
to  high  diffusion  with  low  film  resistance.  At  higher 
values  of  8  and  at  an  overpotential  of  more  than  100 
mV,  the  current  density  approaches  a  limiting  value. 
This  prediction  of  limiting  current  is  possible  only  with 
the  modified  model.  The  previous  model  fails  to  show 
such  limiting  behavior. 

The  effect  of  the  exchange-current  density  (i0)  was 
studied  by  varying  i0  from  0.05  to  1.0  pA  cm-2,  as 
shown  in  Fig.  3.  It  is  found  that  as  i0  increases,  the 
current  generated  by  the  electrode  increases. 

The  grain  radius  ( R )  affects  the  polarization  per¬ 
formance  of  the  electrode,  as  shown  in  Fig.  4  in  which 
R  is  varied  from  15  to  30  pm.  The  larger  the  grain 
radius,  the  less  is  the  current  generated. 

The  effect  of  gas  diffusivity  (D,e)  in  the  electrolyte 
is  shown  in  Fig.  5.  Dle  has  been  varied  from  the  5.0  X 10" 5 
to  20X10"5  cm2  s"1.  As  Dlc  increases,  the  current 
density  increases  for  each  value  of  the  overpotential. 
This  is  expected  because  more  reactant  is  available  for 
the  reaction. 

Fig.  6  shows  the  effect  of  hydrogen  diffusivity  (Dg) 
in  the  macropores  of  the  electrode  for  values  in  the 
range  from  2.5X10”12  to  50X10"12  mol  cm"1  kPa  1 
s'1.  AsDg  increases,  the  current  density  increases.  This 
behaviour  reflects  the  significance  of  the  increase  in 
Dg  in  higher  supply  of  reactant  gas  to  the  electrode. 

The  reactant  pressure  profile  in  the  electrode  is 
shown  in  Fig.  7  for  different  values  of  the  macropore 
diffusivity,  Dg,  in  the  range  from  0.1  X 10  11  to 
1000x10“ 11  mol  cm"1  kPa"1  s"1.  At  higher  values  of 
Dg,  the  pressure  in  the  macropores  remains  almost 
constant,  since  the  macropore  diffusion  resistance  is 
then  negligible  and  the  overall  reactant  is  controlled 
by  the  electrochemical  reaction  and/or  micropore  dif¬ 
fusion.  At  lower  values  of  Dg,  the  overall  reaction  will 
be  macropore-diffusion  controlled.  At  extremely  low 
values,  there  will  be  very  little  reactant  available  for 
the  reaction,  and  a  limiting  current  will  result. 
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Table  1 

Estimated  and  fitted  parameters  for  modified  model  at  different  temperatures 


Temperature 

(°C) 

Estimated  parameters  from  literature 

Parameters  found  from  modified  model 

(n-‘  cm- 

kXlO9 

)  (mol  cm-5  kPa-1) 

DgXl0" 

(mol  cm-  kPa“' 

£>,X  105 

s'1)  (cm2  s“') 

“ 

SxlO5 

(cm) 

/„  x  106 
(mA  cm”2) 

25 

0.254 

1.449 

4.589 

1.515 

0.59 

7.0 

11 

35 

0.300 

1.230 

4.097 

1.938 

0.59 

7.0 

26 

45 

0.344 

1.0145 

3.540 

2.625 

0.59 

7.0 

85 

55 

0.390 

0.8696 

3.182 

3.931 

0.59 

8.0 

160 

65 

0.434 

0.6667 

2.551 

7.404 

0.59 

8.0 

300 

75 

0.464 

0.5070 

2.279 

9.529 

0.59 

8.0 

1000 

Overpotential,  mV 

Fig.  2.  Effect  of  film  thickness  on  electrode  polarization  using 
parameters  at  348  K. 


Overpotential,  mV 


Fig.  4.  Effect  of  grain  radius  on  electrode  polarization  using  parameters 
at  348  K. 


Overpotential,  mV 

Fig.  3.  Effect  of  exchange-current  density  on  electrode  polarization 
using  parameters  at  348  K. 


Overpotential,  mV 

Fig.  5.  Effect  of  diffusivity  of  hydrogen  in  electrolyte  on  electrode 
polarization  using  parameters  at  348  K. 


3.2.  Comparison  of  modified  model  prediction  with 
experimentally  obtained  data 

The  polarization  data  for  hydrogen  oxidation  on  Ni/ 
PTFE  electrodes  in  25%  KOH  electrolyte  at  different 
temperatures  (25  to  75  °C)  has  been  used  for  comparison 


with  data  predicted  by  the  modified  model.  The  pa¬ 
rameters  used  to  solve  the  modified  model  were  obtained 
through  electrode  characterization,  estimation  by  em¬ 
pirical  expressions,  and  published  data.  Conductivity 
of  25%  KOH  solution  was  estimated  by  an  empirical 
formula  [6].  The  Henry’s  law  constant  at  different 
temperatures  was  obtained  from  the  solubility  data 
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Overpotential ,  mV 


Fig.  6.  Effect  of  diffusivity  of  hydrogen  in  macropores  on  electrode 
polarization  using  parameters  at  348  K. 


Fig.  7.  Reactant  pressure  profiles  in  electrode  at  different  macropore 
diffusivities  at  ij0  =  44  mV  and  other  parameters  at  348  K. 

given  by  Shoor  et  al.  [7].  The  diffusivity  of  hydrogen 
in  the  macropores  was  calculated  by  the  Hirschfelder’s 
formula  [8].  The  diffusivity  of  hydrogen  in  KOH  solution 
was  estimated  by  a  relation  given  in  Ref.  [9]  using 
viscosity  and  density  data  given  in  Ref.  [10].  With  these 
values,  the  model  equations  were  solved  and  fitted  to 
the  experimental  data  (the  fitted  intrinsic  parameters, 
namely,  exchange-current  density,  charge-transfer  coef¬ 
ficient  and  film  thickness,  were  obtained  in  this  manner). 
Both  estimated  and  fitted  parameters  are  listed  in  Table 
1.  Besides  these,  the  following  parameters  and  constants 
were  also  used: 

2=1.0 
n  =  2.0 

F= 9.6484531  X 104  C  mol-1 
R  =  2.0X10- 3  cm  (Ref.  [3]) 

L  =  0.032  cm  (Ref.  [3]) 

Rg  =  8.31  J  mol  K”1 


emac  =  0.21  (Ref.  [3]) 
a  =  5.4x10s  cm2  cm-3  (Ref.  [3]) 

Polarization  curves  obtained  from  this  modified 
model,  along  with  the  previous  model  and  experimental 
data  are  plotted  in  Figs.  8  and  9  for  temperatures  of 
35,  45,  55  and  75  °C.  It  is  observed  that  the  modified 
model  fits  the  experimental  data  very  well.  It  can  be 
inferred  from  the  fitted  data  that  the  film  thickness 
remains  constant  with  temperature,  while  the  exchange- 
current  density  increases.  In  fact,  the  exchange-current 
density  is  the  electrochemical  counterpart  of  the  rate 
constant  for  the  chemical  reactions.  Therefore,  an  Ar¬ 
rhenius’  type  of  temperature  dependence  is  expected. 
The  logarithm  of  i0  is  plotted  against  the  inverse  of 
the  absolute  temperature  in  Fig.  10.  The  plot  is  linear 
and  the  estimated  value  of  the  activation  energy  is  75.2 
kJ  mol-1.  This  value  is  about  twice  that  predicted  by 
the  spherical-grain  model  [4].  This  indicates  that  the 
hydrogen  oxidation  on  Ni/PTFE  electrode  is  mass- 
transfer  limited.  Since  the  previous  spherical-grain 
model  did  not  include  mass-transfer  resistance,  the 


Fig.  8.  Comparison  of  experimental  data  with  old  and  modified  model 
(7=308  and  318  K). 


Fig.  9.  Comparison  of  experimental  data  with 
(7=328  and  348  K). 


i  old  and  modified  model 
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Fig.  10.  Exchange-current  density  vs.  1000/7"  (Arrhenius’  plot). 


activation  energy  calculated  was  the  apparent  one.  By 
contrast,  the  modified  model  delivers  the  true  value. 
As  expected  from  mass-transfer-limited  catalytic  re¬ 
actions,  the  apparent  activation  energy  is  about  half 
of  the  true  activation  energy  [11]. 


4.  List  of  symbols 

a  surface  area  of  the  catalyst  (cm2  per  unit  volume) 
A  see  Eqs.  (5)-(7) 

C  concentration  of  reactant  in  the  grain  (mol  cm  “ 3) 

Cc  concentration  of  reactant  in  the  film  (mol  cm-3) 
C0  bulk  concentration  of  reactant  (mol  cm-3) 

Cs  concentration  of  reactant  at  the  interface  of  film 

and  grain  (mol  cm-3) 

D'g  diffusivity  of  reactant  gas  into  the  macropores 
(cm2  s-1) 

Dg  D’g  k  (mol  cm-1  kPa-1  s-1) 

£>,  diffusivity  of  dissolved  gas  into  the  electrolyte 
(cm2  s-1) 

Dle  diffusivity  of  dissolved  gas  into  flooded  grain 
(cm2  s-1) 

F  Faraday’s  constant  (9.6484531  X 104  C  mol-1) 
k  Henry’s  constant 

i  local  current  density  (mA  cm-2) 

t0  exchange-current  density  (mA  cm-2) 

j  measurable  current  density  (mA  cm-2) 

L  thickness  of  the  electrode  (cm) 

n  number  of  electrons 

Ni  number  of  grains  per  unit  surface  area  of  the 
electrode  (cm-2) 
p  pressure  (Pa) 


p0  bulk  gas  pressure  (kPa) 
r  radial  direction  (cm) 

R  grain  radius  (cm) 

R*  consumption  of  reactant  in  the  electrode  (mol 
cm-3  s-1) 

Rg  gas  constant  (8.314510  J  K-1  mol-1) 

rgas  consumption  of  reactants  in  the  grain  (mol  cm  -  3 

T  reaction  temperature  (K) 

x  axial  direction  (cm) 

z  stoichiometric  number 

Greek  letters 

a  charge-transfer  coefficient 

5  film  thickness  (cm) 

7)  overpotential  (mV) 

kc  effective  conductivity  of  the  catalyst  (fl-  1  cm-1) 

emac  macroporosity  of  the  electrode 
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